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DETERMINATION  OF  THE  OPTICAL  CONSTANTS  OF  HIGHLY 
ABSORBING  MATERIALS  IN  THE  INFRARED  REGION 

1.  INTRODUCTION . 

The  description  of  the  interaction  of  light  with  matter, 
as  formulated  from  Maxwell's  electromagnetic  theory  uses 
microscopic  parameters  to  characterize  the  specific  pro¬ 
perties  of  a  material  involved  in  the  process.  These  para¬ 
meters  (notably  the  dielectric  constant)  are  related  to  the 
optical  constants,  the  refractive  index  n  and  the  absorption 
coefficient  k  of  the  material  in  question.  In  a  simplified 
picture,  the  index  n  characterizes  the  refraction  properties 
of  the  material,  since  it  is  the  ratio  of  the  velocity  of 
the  wave  inside  to  that  outside  of  the  material,  and  the  ab¬ 
sorption  coefficient  k  describes  the  attenuation  of  the  wave 
as  it  penetrates  the  medium. 

The  optical  constants  can  be  obtained  by  an  analysis  of 
reflection  and  transmission  data.  One  determines  n  and  k 
using  polarized  light  measurements,  since  the  reflection  and 
transmission  properties  of  the  material  are,  in  general,  dif¬ 
ferent  for  light  polarized  parallel  to  or  perpendicular  to 
the  plane  of  incidence.  For  crystals  with  optical  axes,  the 
optical  constants  depend  upon  the  orientation  of  the  surface 
relative  to  plane  of  incidence.  This  latter  effect  is  lost 
in  polycrystalline  materials  where  only  average  optical  con¬ 
stants  can  be  obtained.  Reducing  the  size  of  the  domains  of 
the  polycrystalline  sample  results,  eventually,  in  an  amor¬ 
phous  material  where  the  long  range  order  of  the  atoms  has 
disappeared  and  only  short  range  order  of  the  molecular 
type  prevails. 

Powders  have  the  nonoriented  properties  of  the  poly¬ 
crystalline  materials  but,  in  addition,  the  powder  particles 
are  now  embedded  in  a  dielectric  which  might  be  air,  KBr, 
other  materials,  or,  in  the  limiting  case,  a  vacuum.  The 
discontinuity  of  the  optical  constants  at  the  surface  of  the 
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powder  particles  introduces  alterations  in  the  transmitted 
or  reflected  light.  One  effect  is  that  light  is  scattered 
out  of  the  direction  of  the  reflected  or  transmitted  light 
beams.  This  scattered  light  is  unavailable  for  analysis; 
its  absence  introduces  an  error  in  the  determination  of 
the  optical  constants.  This  effect  is  particularly  trouble¬ 
some  if  the  absolute  values  of  reflected  and  transmitted  in¬ 
tensities  are  required.  Nevertheless,  using  the  pellet 
method,  in  which  the  powder  is  embedded  in  KBr  and  the  op¬ 
tical  constants  are  obtained  from  transmission  and  reflec¬ 
tion  data1,  surprisingly  good  results  have  been  obtained. 

In  this  work,  the  "lost"  light  is  simply  disregarded  in  de¬ 
termining  the  optical  constants.  Observations  of  the  trans¬ 
mitted  and  reflected  light  from  powder  samples  have  also 
shown,  in  some  instances,  new  spectral  features  when  compared 
to  the  bulk  materials.  Careful  investigation  has  shown  that 
these  features  depend  on  the  size,  shape  and  nature  of  the 
dielectric  in  which  the  powder  particles  are  embedded.  The 
most  interesting  change  occurs  with  the  optically  active 
transverse  mode  frequency  in,  for  example,  alkali  halide 
crystals.  The  frequency  of  this  mode  shifts  to  higher 
values  if  the  size  of  the  crystallite  is  reduced. 2 

Theoretically,  one  can  account  for  these  phenomena 
using  the  phonon  concept.  One  considers  finite  atomic 
chains  but  omits  the  application  of  cyclic  boundary  condi¬ 
tions  used  in  large  crystals.  An  alternate  approach,  the 
continuum  theory,  first  formulated  by  Frohlich3,  considers 

*F.  E.  Volz,  Appl.  Optics  12,  564  (1973) 

2J.  T.  Luxon,  D.  J.  Montgomery  and  R.  Summitt,  Phys.  Rev. 

188,  1345  (1969) 

SH.  Frdhlich,  Theory  of  Dielectrics,  (Oxford  University 

Press,  Oxford,  1948) 
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the  polarizations  in  the  mixed  medium.  The  effect  of  par¬ 
ticle  shape  was  considered  for  spherical  and  cylindrical  geo¬ 
metries.  The  latter  can  be  extended  to  the  case  of  fibers. 
The  general  approach  to  dielectric  constant  studies  also  pro¬ 
vides  the  possibility  of  studying  metal  particles  embedded, 
either  in  their  own  oxide  layers,  or  in  a  dielectric  matrix. 
Before  turning  to  the  more  detailed  discussion  of  powder  sam¬ 
ples  we  will  discuss  the  determination  of  optical  constants 
of  amorphous  bulk  materials. 


2.  SOLUTION  TO  THE  WAVE  EQUATION  FOR  A  LOSSLESS  DIELECTRIC 
We  start  with  Maxwell's  equations  written  as 

SB 


V  x  E  =  - 


V  x  B  =  y 


3 1 
3D 

at 


(1) 

(2) 


D  -  P 


V  •  B  =  0  (4) 

where  E  is  the  electric  field  vector,  D  is  the  displacement 
vector  defined  by  D  =  e  E  where  e  is  the  permittivity,  B  is  the 
magnetic  induction  vector,  y  the  permeability  of  free  space 
and  p  is  the  charge  density  (zero  in  our  case) . 

From  Eqn. (1)  and  Eqn. (2)  we  obtain 

-4- 

VxVxE  =  -VxB  =  - 


(3) 


V  (V-E)  -  72E  =  -  yQ  D 

* 

CE  =  y  e  E 
0 


y  2ir  = 


(5) 


This  is  the  wave  equation  for  the  dielectric  medium  with  per¬ 
mittivity  e.  It  reduces  to  the  vacuum  equation  for  e  -*•  eQ. 

As  solutions  to  the  wave  equation  we  will  consider  only 
plane  waves  propagating  in  the  x  direction 


where  v  is  the  velocity  of  the  wave.  Eqn.(5)  with  E(x,t) 
has  the  form  of  Eqn. (8) .  Consequently,  we  can  take 
yQ  e  =  1/v2  and,  using  the  definition  of  the  refractive 
index,  n  =  c/v  get 

n2  =  c2  e  (9) 

The  permittivity  e  can  be  expressed  as  the  product  of  the 

dielectric  constant  k  and  the  permittivity  of  free  space  eq 

e  —  k  e  (10) 

0 

Introducing  Eqn. (10)  into  Eqn. (9)  gives 

2  2 

nz  =  c4  y  e  k 
o  o 

and,  since  V0e0  =  1/c2  ,  we  have  finally 

n  =  /~ic  (11) 

The  refractive  index,  for  a  non-magnetic  material  (y  =  uQ) , 
is  equal  to  the  square  root  of  the  dielectric  constant. 
Another  useful  relation  is 


(12) 
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3 .  FRESNEL ' S  FORMULAE 

To  obtain  Fresnel's  formulae,  we  consider  the  behavior 
of  the  electric  and  magnetic  field  vectors  at  the  interface 
separating  the  media  with  refractive  indices  n^  and  n2»  The 
boundary  conditions  requires  that  the  tangential  components 
of  both  E  and  B  are  continuous  at  the  interface.  We  will 
adopt  the  notation  of  Fig.  1. 


Wave 

Figure  1.  Fresnel  Formula  Figure 


3.1  E  Perpendicular  to  Plane  of  Incidence 
We  represent  the  incident,  reflected,  and  transmitted 
plane  waves  by 


=  A  e 

(13a) 

E;  =  A 1  e  i(K1x’-&Jt) 

(13b) 

E"  =  A"e  i(K2x"-a)t) 

(13c) 

and  K2  are,  respectively,  the  wave 

vectors  outside 

and  inside  the  dielectric.  Note  that  the  angular  fre¬ 
quency  a)  is  the  same  for  both  media.  From  the  boundary  con¬ 
dition  on  E  we  obtain 

E  +  E'  =  (14) 

Since  the  phase  factors  are  equal  at  the  boundary,  we  have 
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ym 


the  result 

A  +  A'  =  A"  (15) 

The  first  Maxwell  equation  relates  the  E  and  B  vector 
as  follows  - 


7xE  =  -i(Bxi  +  Byj+Bzk)  (16) 

From  the  y  and  z  components  of  the  electric  field 

Ey  =  A  e  (17a) 

E  ■  A  e  l(Kx"ut)  (17b) 

z 

which  are  also  propagating  in  the  x  direction,  we  get, 
using  Eqn. (16) 


3B  3E 

_ z  __ 

3t  "  3x 

Integration  results  in 


3B  3E 


and 


3t  3x 


(18) 


Bz  =  -  iK/Eydt  -  -  iKA 


1KX  je  -i«t  dt 


-iK  Agi  (Kx-ut)  _  K  E  _  n  E 
-  xco  £o  y  c  y 


and  correspondingly 


B  =  -  -  E 
y  c  z 


(19) 


(20) 


where  n  is  either  n^  or  n2* 

The  B  vector  accompanying  the  E,  E'  and  E"  vectors  in 
the  above  figure  are  in  the  plane  of  the  paper  and  have,  at 
the  boundary,  the  components 

-  Bj|  cos  6  ,  B||  cos  9  ’  ,  B"  cos  9  " 

Therefore  the  boundary  condition  on  B  can  be  written  as 


-  B||  cos  9  +  Bj|  cos  9  '  =  -  B”  cos  9"  (21) 
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Using  the  relation  between  the  B  and  E  and  the  fact  that  the 
phase  has  to  be  the  same  for  all  terms  on  the  boundary  we 
get 


-n^  A  cos  0  +  n^  A  '  cos  0‘  =  -  nj  A"  cos 

Eqns. (15)  and  (22)  can  be  written  as 

0  "  (22) 

A 1  A" 

and 

(23) 

1  -  JT  ‘  “  *- 

using 

n_cos  0 " 

(24) 

Cl  E  - 

n^cos  0 

(25) 

We  now  introduce  the  amplitude  reflection 

coefficient 

A 1 

rEj_  "  A 

and  the  amplitude  transmission  coefficient 

(26) 

A" 

= 

Eqns. (23)  and  (24)  can  now  be  written  as 

(27) 

rEj.  '  *ex  "  ‘  1 

(28) 

rEx  +  “V  -  1 

(29) 

which  are  to  be  solved  for  r_  and  t_  .  The  result  is 

Ex  E± 

n1  cos  0  -  n,  cos  0 " 

_  i. _ z _  /  *3  n 

rEj_  n£  cos  0  +  ^  cos  0 "  1 


2n^  cos  0 

n^  cos  0  +  n£  cos  0  " 


(31) 
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3.2  E  Parallel  to  Plane  of  Incidence 


In  a  similar  way,  the  amplitude  reflection  coefficient 
r^  and  amplitude  transmission  coefficient  t^  for  the  case 
where  the  E  is  in  the  plane  of  incidence,  are  given  by 


n^cos  8"  -  njcos  0 
n^cos  8*'  +  n2cos  0 


2n^cos  0 

n^cos  8"  +  n2cos  0 


4.  DETERMINATION  OF  REFRACTIVE  INDEX 

We  consider  an  air-dielectric  interface  and  take  n^  =  1 
and  n2  =  n.  At  normal  incidence  we  have  the  amplitude  ratios. 

r  >  (34) 


*  *  TTF  1351 

However,  experimentally,  we  measure,  for  example,  the  ratio 
reflected  power  flux  „  _ 


incident  power  flux 


=  R  =  Reflectance 


Using  Poynting's  vector,  the  reflected  power  flux  is  given 
by 

R  =  r2  -  (~)  2  (36) 


Similarly ,  one  obtains  for  the  transmitted  power  flux 

n  n  x  n 


T  =  at2  = 


2  cos  9" 
i,  cos  0 


A »  _ 

(•V)2 


At  normal  incidence,  we  can  obtain  the  refractive  index 
from  one  measurement  of  R,  that  is  Eqns.(34)  and  (36)  can  be 
combined  to  give 


1  -  r  1 


For  non-normal  incidence  we  can  write  the  reflection  co¬ 


efficients  as 


r  =  COS  9  -  n/1  -  (sin  9/n)z 
EX  cos  9  +  n/1  -  (sin  0/n)2 

and 

-n  cos  0  +  /I  -  (sin  0/n)2 

rE  =  - ; - - - —  (40) 

||  n  cos  +  /l  -  (sin  0/n)2 

using  the  law  of  refraction.  One  sees  again  that  a  single 
measurement  is  sufficient  to  determine  n  if  the  angle  of  in 
cidence  is  known  and  polarized  light  is  used. 


5.  DISPERSION  MODEL  OF  A  DIELECTRIC 

In  another  report  we  investigated  the  relationship  be 
tween  the  refractive  index  and  the  parameters  of  an  oscil¬ 
lator  model  having  resonance  frequencies  (do^,  mass  nu, 
charge  e  and  damping  constant  y.  If  the  lossless  case  is 
considered,  that  is  no  damping,  the  permittivity  is  found 
to  be 


^-  =  (n*)2  =1  +  1 


eomi(a3oi-<0  > 


- 


•2  2 
1  U)  GO* 


where  (o^is  the  plasma  frequency  given  by  oipJ2  =  NPe2/m^eo. 
For  a  real  dielectric,  there  is  always  a  polarization  ef¬ 
fect  from  the  surrounding  medium.  At  a  particular  point, 
the  so  called  Lorentz  contribution  to  the  polarization  must 
be  taken  into  account.  This  results  in  a  shift  of  the  re¬ 
sonance  frequencies  to  new  values  represented  by  oPT  . 

Let  us  consider  a  system  having  only  two  resonance 
frequencies  and  co^  .  Take  <o  ^  in  the  infrared  region 
and  <nQ 2  at  a  much  higher  frequency.  We  then  have 


=  1  + 


- 2  2 

0)  ~  -  (ti 

o2 


15 


The  last  term  in  the  above  equation  represents  the  high  fre¬ 
quency  electronic  polarization.  Its  contribution  to  e  in 

the  infrared  is  constant.  If  we  define  e(®)  =  e  we  can 

00 

write  „  ...  2 


£  (°°)  _ 
e _  ~  e. 


0> 


=  1  + 


£2_ 


0  "0 

then  Eqn. (42)  becomes 


- 2  2 

“o2  '  “ 


£  (01) 
e 


01 


EL 


ol 


2  2 
-  CD 


(43) 


(44) 


We  can  also  introduce  the  static  permittivity  e(o)=e  and 


get 


£  ( o )  -  e„ 


0)  ,  ‘ 

-  . Rk 


(45) 


6) 


Ol 


Substituting  “pl2  from  this  expression  into  Eqn. (44)  gives 


0) 


e(w)  =  £„  +  (e0-£j 


ol 


(46) 


From  this  expression  we  see  that  if  we  go  from  the 
infrared  resonant  frequency  to  e  (cu) 
lower  frequencies,  the  per¬ 
mittivity  approaches  the 
static  value.  If  we  go  to 
higher  frequencies  the  per¬ 
mittivity  approaches  the 
electronic  value.  The  re¬ 
fractive  index  and  permit¬ 
tivity  are  connected  by  n2  =  e. 

Therefore,  the  physical  origin 
of  n  at  frequencies  above  and 
below  the  infrared  region  can 
be  attributed  to  these  two  factors. 


Figure  2.  Permittivity  and 
Resonance  Frequency 
Relationships 
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6.  STRONGLY  ABSORBING  DIELECTRIC  MEDIA 

We  now  want  to  consider  strongly  absorbing  media,  which 
are  of  prime  interest.  We  begin  with  Maxwell's  equations 
for  the  general  case  of  dielectric  and  conductive  medium. 


V  x  E  =  - 


7  x  B  = 


V  •  $  - 


e  y  If-  +  y  If  +  y  J 

o  o  3t  o  3t  o 


V  •  B  =  0  (50) 

In  addition  to  the  previously  defined  quantities,  we  now 
have  the  polarization  vector  P  =  x*£QE/  where  x*  is  the  com¬ 
plex  electric  susceptability ,  and  the  current  density  vec¬ 
tor  J  =  oE,  where  a  is  the  conductivity.  A  manipulation  of 
Maxwell's  equations,  very  similar  to  that  done  for  the  loss¬ 
less  dielectric  case,  results  in  the  following  wave 
equation 

V’g  =  y  e  221+  y  y  —  (51) 

0  °  at2  0  at*  0  at 


If  the  dielectric  properties  of  the  material  are  re¬ 
presented  by  one-dimensional  oscillator  model  with  losses  , 


we  can  write  P  =  e  x*E  =  e  cu  zE/(<o  2-wz-iyo>) 

o  p  o 


If  the  con¬ 


ductive  properties  of  the  medium  are  represented  by  an  os¬ 
cillator  model  without  restoring  forces,  we  can  write 


J  =  -  E  =  a*E 

1-1(0T 

The  wave  equation  may  then  be  written  as 


y  e  -  +  y  c  x 

0  0  8tz  0  0 


*  9  E  .  _*3E 

*  —  +  y  0* - 

^  J-  2  0  rv  . 
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r 


Introduction  of  a  trial  solution  E  *  E  e  *  ^Kx  "  into 

o 

Eqn. (53)  gives 


or 


K2  *  y  e  to2  +  y  e  X*^2  ~  itoy  a* 
o  o  o  oA  o 


(54) 


K  2 

—  =  y0EQ  +  yoeox*  -  iyQc*/(o  -  yQe* 

Comparing  Eqn. (54)  with  its  lossless  analogue,  Eqn. (12) ,  we 
see  that  for  this  case,  we  have  an  effective  complex  per¬ 
mittivity  e*  instead  of  the  real  permittivity  e.  In  Eqn. 
(12) ,  for  the  lossless  case  we  saw  that  the  refractive  in¬ 
dex  and  permittivity  are  related  by  n  =  /e/e  .  A  complex 
e  implies  a  complex  n,  thus  n*  =  /e*/eQ  now  describes  the 
optical  properties  of  the  material.  If  we  represent  the 
complex  refractive  index  by  n*  =  n  +  ik,  we  see  that  the 
losses  are  associated  with  the  imaginary  part  k,  called  the 
extinction  coefficient. 

The  formal  substitution  of  e*  for  e  and  n*  for  n  does 
not  change  the  formal  derivation  of  the  law  of  refraction 
and  Fresnel's  formulae,  the  law  of  refraction  becomes 


sin  6  =  n*  sin  0"  (55) 

and  the  reflection  coefficients,  obtained  from  Fresnel's 
formulae  are,  for  n^  =  1, 


and 


cos  9  -  n2*  cos  9" 

£  — 

X  cos  0  +  n^*  cos  9" 

cos  0”  -  n,*  cos  0 

'E  - - - 

II  cos  9"  +  n2*  cos  0 


(56) 


(57) 


There  are  three  points  to  be  mentioned  here: 

1)  Both  r„  and  r_  are  complex  and  can  be  written 

Ei.  Eti 

as  a  product  of  an  amplitude  factor  and  a  phase 
factor  as 
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r  =  p  e  iaJ.  and  rp  =  p  e  ia\\ 
h±  ±  *  II 

The  angle  a  is  the  phase  shift  introduced  after 
reflection. 

2)  There  is  no  Brewster  angle  for  r„  .  We  can  see 

eii 

this  by  introducing  n*2  =  n2  + ik2  into  the  Brewster 
angle  relation  tan6fi  =  n2*  .  Then  tane0  =  n2  +  ik2 . 
Since  e0  is  real,  no  solution  exists. 

3)  There  is  an  attenuated  wave  in  the  absorbing 
medium.  Considering  only  normal  incidence, 
we  have 

E  =  E  e  i <***-“*>  =  E  e  il|(n+ik)x-03t] 

o  o 

E  =  (e  "  £kx)  Eoe  i(Hnx"“t) 

where  the  first  factor  on  the  right  side  is  the 
attenuation  factor  and  the  second  represents  a 
travelling  wave.  The  amplitude  of  the  wave  is 
attenuated  to  a  value  of  1/e  after  travelling  a 
distance 

d-  =-!_=_£_=  _ *  _L_ 

u  „  03  K  2ircK/X  2ttK 

c  K 

into  the  medium.  Since  the  intensity  is  propor¬ 
tional  to  | E 2 | ,  the  power  flux  decreases  to  1/e  of 
its  value  at  the  surface  in  a  distance  d  =  X/4iTk. 

We  see  that  the  optical  constant  k  can  be  deter¬ 
mined  through  attenuation  measurements. 

6.1  Evaluation  of  Optical  Constants  of  Highly 
Absorbing  Materials.  Graphical  Methods 
The  reflection  coefficients,  as  given  in  Eqns. (56)  and 
(57) ,  can  be  expressed  as  functions  of  the  complex  refractive 
index  n*  and  the  angle  of  incidence  6  if  Snell's  law  Eqn. 
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(55),  is  used  to  eliminate  9".  If  we  calculate  the  square 
of  the  absolute  value  of  the  reflection  coefficients,  we 
get 

rE  I’  =  I.  =  (a  -  cos  b) 2  +  b*  ,58, 

EJ-  (a  +  cosb)2+b2 

and 


2  _ 


=  I. 


II 


where 


(a  -  sin  9  tan  9) 2  +  b2 
(a  +  sin  9  tan  9)  2  +  b2 


(59) 


a2  ■  [ (n2-k2-sin29) 2+4n2k2]^+(n2-k2-sin29) } 


b2  =j{  [  (n2-k2-sin29)  2+4n2k2]i- (n2-k2-sin29) } 


These  complicated  formulae  relate  the  intensity  ratios  of 
reflected  light,  for  parallel  and  perpendicular  components, 
to  the  angle  of  incidence  9  and  the  optical  constants.  To 
determine  the  optical  constants,  one  needs  the  inverse  re¬ 
lationship,  that  is,  n  and  k  as  a  function  of  the  ratios  of 
the  reflected  light  and  the  angle  of  incidence.  Unfortu¬ 
nately,  the  above  equations  can  not  be  inverted  analyti¬ 
cally.  However,  J.  Vasicek  has  calculated  "Tables  of 
Determination  of  Optical  Constants  from  the  Intensities  of 
Reflected  Light",  (Czechoslovakian  Academy  of  Science, 
Prague,  1963) .  In  these  tables,  the  following  measurable 
quantities  are  tabulated  as  a  function  of  9,  n  and  k: 


Numerical 


values  are  given  for  0  i  9  $  80°,  in  increments  of  5°. 
Values  of  n  and  k  are  given  in  increments  of  0.2  for 
0.2  5  n  £  7.0  and  0.0  5  k  £  6.0.  A  graphical  analysis  of 
the  reflection  data  using  Vasecek's  tables  can  be  used  to 
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^sample 

Figure  3.  Intersection  c*  the 
Three  Isoreflectance  Curves 


give  values  of  n  and  k  as  follows:  For  a  given  angle  of 

incidence  one  plots  on  an  n-k  curve,  the  locus  of  all  points 

having  the  same  measured  re-  9 

fleeted  intensity  I_  ,  I_  ,  n 

e-L  e||  - _  yT 

I  or  Ip.  This  procedure  n 

is  repeated  for  a  second,  sample  \ 

and  preferably  a  third,  y'  i  ^ 

angle  of  incidence.  The  '  I 

intersection  of  the  three  1 

i 

isoreflectance  curves  gives  i 

the  optical  constants  of  the  , 

material,  see  Fig.  3.  ^sample 

Normally  the  curves  do  not  Figure  3.  Intersection  o*: 

intersect  at  a  single  point,  Three  I=°re£lectance  Curv 

but  define  a  triangularly 

shaped  region  which,  in  effect,  is  an  indication  of  the 
inherent  error  in  the  method.  For  such  cases,  the  'center 
of  gravity'  of  the  triangle  is  used  to  extract  the  values 
of  n  and  k  for  the  sample. 

The  behavior  of  tne  reflected  intensities,  as  a  func¬ 
tion  of  angle  of  incidence,  can  be  seen  by  plotting  curves 
of  I„  and  I_  as  functions  of  0.  The  parameters  affect- 
ing  the  shapes  of  the  curves  are  n  and  k.  W.  R.  Hunter1* 

has  generated  such  plots  for  0.2  1  k  -  3.3  in  increments  of 

0.1  and  0.3  S  n  $  2.3  in  increments  of  0.5.  A  typical  plot 

for  n=1.3  and  k  =  2.3  is  shown  in  Fig.  4. 

The  difference  between  Fig.  4  and  the  corresponding 
plot  for  a  lossless  (k  *  0)  dielectric  is  that  for  the 
lossless  case,  I_  =  0  at  the  Brewster  angle.  When 


k  >  0,  I, 


never  goes  to  zero. 


*  W.  R.  Hunter,  J.  Opt.  Soc.  Am.  55  1197  (1965) 
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In  Fig.  4,  we  see 


that  I. 


and  I. 


“II  "J. 

differ  most  for  large 
values  of  8.  Thus, 
experimentally,  it  is 
desirable  to  choose 
appropriately  large 
values  of  0  to  obtain 


Figure  4 .  1  Values  at 

Different  Angles 


the  most  accurate  values  of  n  and  k. 

6.2  Ellipsometry 

The  reflection  measurements  considered  so  far  have  re¬ 
quired  the  use  of  absolute  values  of  two  reflected  inten¬ 
sities  or  the  ratio  of  two  measured  reflection  values.  If 
light  is  scattered  out  of  the  light  path  and  thus  not  ac¬ 
counted  in  the  determination  of  the  reflected  intensities, 
errors  in  the  values  of  n  and  k  will  result.  Since  the  pow¬ 
der  particles  have  a  random  orientation  with  respect  to  the 
plane  of  incidence,  one  may  assume,  in  the  first  approxi¬ 
mation,  that  the  lost  scattered  light  is  the  same  for  both 
components  of  the  reflected  light.  If  one  derives  formulae 
from  the  original  Fresnel's  equations,  which  make  use  of  the 
difference  between  the  intensities  of  the  polarized  compo¬ 
nents,  the  error  cancels  out.  If  a  ratio  of  intensities  is 
used,  the  error  is  reduced. 

Ellipsometry,  which  has  been  used  extensively  in  the 
visible  region  to  determine  optical  constants  of  metals, 
uses  only  one  ratio  of  reflection  measurements.  The  second 
required  piece  of  information  follows  from  the  fact  the 


phase  shift  differences  between  r^  and  rr 
principal  angle,  is  90 
polarized. 


at  the 


‘1 

That  is,  the  light  is  circularly 


If  one  determines  the  principal  angle  9^  and 
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then  measures  the  ratio  | r_  /r  |  at  6  ,  one  will  have 

E||  P 

the  necessary  experimental  information  to  calculate  n  and  k. 
The  principal  angle  is  determined  using  a  compensator. 

Using  the  phase  difference  property  of  r„  and  r„  ,  Eqns. 

EX 

(56)  and  (57)  can  be  divided  to  yield 


(rE  i  \  cos  6  -  n  *  /l  -  (sin  0/n  *) 2 

- -  = . . . —  — - - -  (60) 

Ell  /  0=6p  /l  -  (sin  0/n2*)2  -  n2*  cos  0 

If  we  calculate  the  square  of  the  absolute  value  of  Eqn. (60) 

we  can,  after  lengthy  calculations,  obtain  the  following  use 

ful  relations  between  the  optical  constants  n  and  k  and  the 

experimental  quantities  0  and  ip  =  I  (r_,  /r_  )D  I2 

Ex  Ejj  0p 

n2  +  k2  =  sin20  +  cos24ip  tan20  sin20 

2nk  =  sin2  4 ip  tan20  sin20  (61) 

Using  a  C02  laser  operating  at  10.6ym,  ellipsometry  has 
been  used  to  evaluate  infrared  optical  constants.  The  un¬ 
availability  of  suitably  transparent  compensator  materials 
over  a  broad  frequency  band  causes  some  experimental  dif¬ 
ficulties.  For  example,  sapphire  can  only  be  used  in  the 
2-3um  region.  Another  consideration  in  the  infrared  is, 
the  determination  of  the  principal  angle.  Using  a  laser, 
with  its  narrow  beam,  no  problem  exists.  But  for  a  non¬ 
laser  source,  to  get  sufficient  energy,  the  beam  diameter 
must  be  enlarged.  The  result  is  a  loss  of  accuracy  in  the 
principal  angle  measurement.  This,  in  turn,  affects  the 
accuracy  of  the  value  of 

We  will  not  discuss  further  the  possible  application  of 
ellipsometry  to  our  problem.  Instead  we  will  concentrate  on 
the  applications  using  (1)  reflection  measurements  at  dif¬ 
ferent  angles,  (2)  reflection  measurements  with  polarized 
light  at  a  single  angle  or  combination  thereof. 


7.  EFFECT  OF  POLARIZATION  PROPERTIES  ON  MEASURED 
REFLECTANCES 

We  assume  that  the  experimental  apparatus  consists  of 
a  conventional  double  beam  infrared  spectrophotometer  with 
an  auxiliary  reflection  unit  mounted  in  the  sample  beam 
light  path.  Of  interest  here  is  the  influence  of  the  po¬ 
larization  properties  of  the  instrument  and  their  affect 
on  the  resulting  measurements.  Initially,  we  will  consider 
only  the  use  of  a  single  reflection  attachment. 

Since  both  reflection  coefficients  are  functions  of  0, 
n,  and  k,  in  general,  measurements  at  two  angles  of  inci¬ 
dence  will  be  required.  That  is,  since  r„  =  f(0,n,k)  and 
r  =  g  (6  ,n,k)  ,  then  f  (6^^  ,n,k) ,  f{02,n,k),  g(91,n,k)  and 
g(02,n,k)  are  required  to  extract  n  and  k  values  from  the 
measured  data. 

A  schematic  diagram  of  the  optical  path  of  the  instru¬ 
ment,  set  up  for  reflection  measurements,  is  shown  in  Fig.  5. 

7.1  Unpolarized  Light.  Two  Angles  of  Incidence 

Let  us  call  the  intensity  of  light  emitted  from  the 
source,  and  polarized  perpendicular  and  parallel  to  the 
plane  of  incidence,  A  and  B  respectively.  Then  A  +  B  re¬ 
presents  the  intensity  of  the  sample  beam.  We  call  A '  +  B ' 
the  corresponding  light  passing  through  the  reference  beam. 
From  the  sample  emerges  light  attenuated  by  a  factor  of  RA 
for  the  perpendicular  component  and  Rn  for  the  parallel  com- 

D 

ponent.  At  the  point  of  recombination  of  the  two  beams, 
the  total  intensity  is  A  +  B  Rg  .  This  amount  of  light 
is  compared  with  the  light  passing  through  the  reference 
beam  A'  +  B'  in  such  a  way  that  a  comb  attenuates  the  re¬ 
ference  beam  until  both  beams  are  of  equal  intensity.  That 
is,  until  we  have 

A  RA  +  B  RA  =  (A*  +  B’)  11  (62) 
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Monochromator/Detector  Sample  Area  Source  Area 


where  presents  the  deflection  of  the  pen  on  the  chart 
paper.  In  order  to  obtain  the  attenuation  of  the  two  po¬ 
larized  beams,  that  is  and  Rg  ,  we  make  the  same  meas¬ 
urement  with  a  standard  replacing  the  sample.  If  it  is 
assumed  that  the  reflectivity  of  the  standard  is  essentially 
unity,  we  obtain  for  that  measurement 

A  +  B  =  (A'  +  B')  T2  (63) 

where  T2  is  the  chart  pen  deflection  when  the  standard  is 
used.  The  desired  result,  that  is  the  reflection  proper¬ 
ties  of  the  sample  are  obtained  from  the  ratio. 


T  AR  +  BRq 

— £_  =  - £ - E_  (64) 

t2  A  +  B  v  1 

The  reflection  coefficient,  using  unpolarized  light  can  only 
be  extracted  if  we  assume  that  A  =  B.  Then  Eqn.  (64)  gives 


RA  +  *B 


(65) 


The  assumption  that  A  =  B  may  be  made  in  most  cases  since 
the  source  emits  unpolarized  radiation  and,  to  a  first 
approximation,  it  may  be  assumed  that  mirrors  do  not  po¬ 
larize  infrared  light. 

We  now  must  account  for  the  polarization  properties  of 
the  instrument  between  the  beam  recombination  point  and  the 
detector.  We  know,  for  example,  that  the  diffraction  grat¬ 
ing  polarizes  the  incident  parallel  and  perpendicular  beams 
differently.  Light  from  the  sample  beam  which  reaches  the 
detector  may  be  represented  by  (ARARA  +  BRgRg)  where  RA  and 
Rg  are  the  attenuation  factors,  due  to  the  instrument  po¬ 
larization  for  the  perpendicular  and  parallel  components, 
respectively.  For  the  reference  beam,  the  corresponding 
term  is  (A'RA  +  B' 

surface,  we  would  get  (ARft  +  BRg)  for  the  sample  beam  and 


Rg) .  For  the  measurement  on  the  standard 
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A 

(A'Ra  +  B'Rg)  for  the  reference  beam.  When  the  beams  are 
combined,  the  equalization  is  given  by 

(ARaRa  +  BRgRg)  =  (A’Ra  +  B'Rb)T1  (66) 

and 

(ARa  +  BRg)  =  A'Ra  +  B'Rq)T2  (67) 

Dividing  Eqn. (66)  by  (67)  gives 


T1  _  arara  +  bVb 

T  A,  /V 

2  ARa  +  BRg 

To  get  the  desired  result,  T1/T2  =  ra  +  RB  ' 


(68) 

we  must  not  only 


assume  that  A  =  B,  as  before,  but  also  that  RA  =  Rg .  That  is, 
the  instrument  polarization  for  each  component  must  be  as¬ 
sumed  to  be  equal. 

Notice  that  both  Eqns. (65)  and  the  result  from  (68) , 
contain  two  unknown  quantities  RA  and  Rg.  To  extract 
either  reflection  component  requires  two  measurements  of 
(Ti/T2)  at  different  angles  of  incidence. 

7.2  One  Component  of  Polarized  Light.  Two  Angles 
of  Incidence 

If  a  properly  oriented  polarizer  is  inserted  into  the 
optical  beam  before  the  detector,  only  one  component  of  the 
total  recombined  beam  will  be  sensed.  If  we  retrace  the 
analysis  of  the  previous  section,  it  is  clear  that  this  case 
can  be  treated  by  setting  either  A  =  0  or  B  =  0.  For  ex¬ 
ample,  if  B  =  0,  Eqn. (68)  gives 


AR.  R. 
A  A 

A 

AR. 


=  R. 


(69) 


Thus,  we  get  the  attenuation  of  one  component  of  the  light 
irrespective  of  the  polarization  properties  of  the  instru¬ 
ment.  Because  the  ratio  T^/T2  is  obtained  by  a  manual  divi¬ 
sion  as  opposed  to  the  ratio  recording  done  by  the  instru¬ 
ment,  possible  differences  between  A  and  A*  are  divided  out 
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and  have  no  effect  on  the  results.  Such  measurements  at 
two  angles  of  incidence  will  yield  the  information  re¬ 
quired  to  determine  n  and  k. 

7 . 3  Both  Components  of  Polarized  Light.  One  Ang le  of 
incidence 

For  each  component  of  polarized  light,  the  argument  pre 
sented  above  in  paragraph  7.2  is  valid.  Since  two  compo¬ 
nents  are  measured,  only  one  angle  of  incidence  is  necessary 

7.4  Ratio  of  Polarized  Liaht  Components.  Two  Angles 


of  Incidence;  No  Standard 


If  measurements  of  the  attenuation  of  reflected  light 
from  the  sample  are  made  for  both  components  of  polarized 
radiation,  and  no  standard  is  used  for  comparison,  the 
equalization  of  the  beam  intensity  is  represented  by 

AR„ R.  =  A ' R. T.  (70) 

A  A  A  A 


bVb  ■  B,RBTB  1711 

Here  is  the  chart  pen  deflection  for  the  sample  with  com¬ 
ponent  A  of  polarized  light  and  B  that  of  component  B.  Note 
that  these  are  two  distinct  measurements  and  that  no  stan¬ 
dard  is  used  for  reference.  The  ratio  Ta/T0  is  given  by 

Tn  AR.B' 

"V  “  brba’  (72) 

and  only  for  A  =  A'  and  B  =  B*  is  the  desired  result  T^/T^  = 


Vrb  • 


obtained.  If  this  measurement  is  made  at  two 


angles  of  incidence,  n  and  k  can  be  extracted  from  the  data. 

Since  we  are  considering  the  use  of  a  double  beam  in¬ 
strument  in  which  the  ratio  recording  method  senses  the  dif¬ 
ference  between  the  sample  beam  and  reference  beam,  it  seems 
reasonable  to  consider  using  a  reflection  unit  in  each  beam. 
In  the  following,  we  consider  arguments  similar  to  those 
presented  above,  but  now  for  the  case  where  the  instrument 
is  operated  with  two  reflection  units. 


7.5  Sample  in  Sample  Beam,  Reference  in  Reference 
Beam.  Unpolarized  light 

For  this  case,  the  equalization  of  the  beam  intensity 
at  the  comb  is  given  by 

arA  +  bVb  *  +  b'*b>  t  (73) 

The  attenuation  due  to  the  sample,  i.e.  T  =  Rft  +  Rg  ,  is 
given  only  if  A  =  A*  =  B  =  B'  and  R.  =  Rn  .  In  addition,  it  is 
assumed  that  the  light  paths  through  both  reflection  units 
are  identical. 

7.6  Sample  in  Sample  Beam,  Reference  in  Reference 
Beam.  One  Component  of  Polarized  Light 

With  a  polarizer  mounted  in  the  common  beam  before  the 

detector,  the  attenuation  of  one  component  (e.g.  B  = 0)  of 

reflected  light  is  given  by 

AR.  R,.  ■  A '  Rn  T  (74) 

A  A  A 

One  only  need  require  that  A  =  A '  to  get  T  =  RA  .  Thus  assum¬ 
ing  equivalent  optical  paths  for  each  beam,  we  can  directly 
get  the  reflection  properties.  Such  measurements  at  two 
angles  of  incidence  will  allow  for  the  calculation  of  n 
and  k. 

7. 7  Samples  in  Both  Beams.  Polarizers  in  Both  Beams 

\  Consider  the  possibility  of  using  identical  reflection 

units,  polarizers  and  samples  in  each  beam.  Then  if  the 
polarizers  are  crossed,  i.e.  oriented  orthogonally,  the 
condition  for  beam  intensity  equalization  is  given  by 

abA  -  B'W  1751 

That  is,  the  polarizer  orientations  produce  A*  =  B  =  0.  The 
desired  result,  T  =  R./R_  ,  is  only  obtained  if  A  =  B '  and 

/V  y\  A  O 

Ra  =  Rg  .  These  restrictions  are  more  seveire  than  those  pre¬ 
viously  discussed.  While  appearing  to  be  an  efficient 
technique,  this  method  actually  is  the  most  difficult  to 
apply. 
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From  the  above  discussions  one  can  see  that  using 
polarized  light  for  measuring  the  reflected  intensities 
from  both  sample  and  standard  require  the  minimum  number  of 
assumptions  on  the  instrument  properties.  For  these  cases, 
it  is  not  necessary  to  assume  equality  of  intensity  of 
source  components  or  grating  polarization. 

From  a  practical  point  of  view,  the  most  convenient 
method  for  isotropic  samples,  for  example  pressed  powder 
disks,  is  to  measure  the  reflected  intensity  from  both 
sample  and  standard  at  one  angle  of  incidence  and  one  di¬ 
rection  of  polarization.  The  measurements  are  then  re¬ 
peated  for  the  other  direction  of  polarization. 

8.  ANISOTROPIC  CRYSTALS 

So  far  we  have  been  concerned  with  the  characteristics 
of  isotropic  materials.  For  anisotropic  crystals,  addi¬ 
tional  considerations  are  introduced.  Recall  that  Fresnel's 
formulae  are  usually  derived  for  isotropic  materials  where 
the  optical  properties  are  the  same  in  all  directions.  For 
uniaxial  crystals,  the  index  of  refraction  has  a  different 
value  along  the  optical  axis  than  it  has  for  the  two  di¬ 
rections  perpendicular  to  the  optic  axis.  To  understand 
this,  we  might  visualize  a  crystal  with  a  structure  having 
planar  symmetry.  That  is,  the  crystal’s  molecular  layers 
are  shaped  such  that  there  is  no  asymmentry  for  the  two 
perpendicular  directions  in  the  plane  of  the  molecules 
(x,y  planes) .  If  the  molecular  layers  are  now  stacked  up, 
as  for  example  with  pancakes,  the  molecular  structure  in  the 
direction  perpendicular  to  the  planes  (z  axis)  is  different 
from  that  in  the  planes.  If  an  electromagnetic  wave  pro¬ 
pagates  along  the  z  axis,  both  electromagnetic  field  com¬ 
ponents  oscillate  in  the  plane  of  the  molecule,  and,  since 
there  is  molecular  symmetry,  they  "feel"  the  same  polari¬ 
zation.  Their  propagation  velocity  is  the  same  and  so  also 


is  their  dielectric  constant  or  index  of  refraction.  For  a 
wave  traveling  along  the  x  or  y  axis,  the  situation  is  dif¬ 
ferent.  The  electric  field  component  parallel  to  the  Z  axis 
will  "feel"  a  different  polarization  of  the  molecules  than 
the  component  in  the  plane  of  the  molecules.  Consequently, 
the  velocity  of  propagation  will  be  different  and  the  in¬ 
dex  of  refraction  will  be  different. 

8 . 1  Measurements  on  Anisotropic  Crystals 
Uniaxial  crystals  are  examined  experimentally  as  fol¬ 
lows.  One  prepares  a  sample  having  its  optic  axis  in  the 
plane  of  the  sample.  For  the  reflection  measurements,  one 
orients  the  sample  in  such  a  way  that  the  optic  axis  is 
oriented  perpendicular  to  the  plane  of  incidence.  In  order 
to  measure  the  optical  constants  parallel  to  the  crystal 
axis,  one  polarizes  the  incident  light  perpendicular  to  the 
plane  of  incidence.  Then  the  sample  is  rotated  by  90  degrees 
and  the  measurements  are  repeated.  For  this  second  case, 
the  incident  light  is  polarized  perpendicular  to  the  optic 
axis.  By  using  polarized  light  on  an  oriented  crystal,  one 
can  study  separately  the  optical  properties  of  the  crystal 
in  directions  parallel  and  perpendicular  to  the  optic  axis. 

9.  APPLICATIONS  TO  POWDER  SPECTROSCOPY 

From  the  above  considerations,  it  is  clear  that  meas¬ 
urement  of  the  optical  constants  of  anisotropic  crystals 
must  be  made  with  care.  The  relationship  between  the  optic 
axis  and  direction  of  incident  light  is  critical.  Conse¬ 
quently,  the  sample  must  be  prepared  appropriately  if  ac¬ 
curate  material  parameters  are  to  be  obtained.  A  similar 
situation  is  found  for  powders.  The  determination  of  the 
optical  constants  of  powders  requires  careful  sample  pre¬ 
paration  and  controlled  measurements.  We  first  discuss 
powders  made  from  isotropic  materials.  It  is  necessary  to 
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consider  refraction  and  scattering  effects.  If  trans¬ 
mission  measurements  are  to  be  made,  the  initial  problem  in 
preparing  a  sample  is  to  eliminate  the  voids  between  the 
powder  particles.  Such  voids  will  cause  the  incident  radi¬ 
ation  to  interact  differently  at  particle-air  interfaces. 
These  different  interactions  will  result  in  an  optical 
characterization  which  can  be  significantly  different  from 
that  of  the  bulk  material.  One  can  try  to  form  a  pressed 
disk  sample  by  subjecting  the  powder  to  a  high  pressure, 
say  in  a  die.  If  the  technique  is  successful  and  the 
sample's  mechanical  properties  are  accomodating,  pressed 
disks  having  a  density  approaching  that  of  the  bulk  can  be 
prepared.  In  such  cases,  the  measured  optical  properties 
should  be  quite  close  to  the  values  determined  from  a  solid. 

An  example  of  such  an  isotropic  material  is  KBr. 

Pressed  KBr  pellets  having  the  visual  transparency  of  a 
polished  crystal  can  be  prepared,  although  probably  not  by 
everyone  and  not  consistently. 

Another  isotropic  material  of  interest  is  ZnS.  East¬ 
man  Kodak  manufactures  optical  quality  disks  obtained  by 
pressing  ZnS  under  high  pressure  and  elevated  temperatures. 
Irtran  2,  as  the  product  is  known  commercially,  is  trans¬ 
parent  to  about  14.5pm.  Its  refractive  index  at  1.5ym 
compares  very  well  to  that  obtained  from  amorphous  films 
and  the  bulk  as  the  following  table  shows. 
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Index  of  Refraction  of  ZnS  at  1.5  ym 


n 

Reference 

Irtran  2 

2.27 

Kodak  Publication  U-72 

amorphous  ZnS  film 

2.26 

J.F.Hall  &  W.F.C.Fergerson 
J.Opt . Soc.Am.  45,714(1955) 

bulk  ZnS 

2.28 

J.R.DeVore, J.Opt.Soc. 

Am.  41,416(1951) 

Well  characterized  materials  such  as  ZnS,  which  are  avail¬ 
able  in  various  forms,  are  of  use  in  studying  the  optical 
properties  of  powders.  Principally,  they  serve  as  "opti¬ 
cal  standards"  by  which  the  accuracy  and  sensitivity 
of  a  given  experimental  method  can  be  checked. 

9.1  Highly  Absorbing  Materials 

For  highly  absorbing  materials,  of  prime  interest  in 
this  work,  transmission  measurements  may  not  be  feasible. 

The  low  concentrations  and/or  thin  sample  paths  required 
for  adequate  energy  throughput  may  introduce  errors  due  to 
unabsorbed  source  radiation  reaching  the  detector.  In  such 
cases,  reflection  measurements  on  solid  samples  must  be  con¬ 
sidered.  Here  again,  surface  conditions  are  important.  For 
sensitive  specular  reflection  measurements,  the  reflecting 
surface  should  not  introduce  significant  scattering  ef¬ 
fects.  We  have  found  that,  for  example,  carbon  and  gra¬ 
phite  powders  can  both  be  cold  pressed  at  room  temperature 
in  a  polished  stainless  steel  die  to  yield  acceptable  re¬ 
flecting  surfaces.  The  graphite  disks  are  quite  stable 
mechanically  while  the  carbon  disks  are  fragile  and  must  be 
handled  with  care.  Providing  the  mechanical  stability  of 
these  pressed  disks  can  be  adequately  preserved,  reflection 
measurements  from  good  quality  surfaces  should  be  realized 
and  the  optical  constants  calculated  from  procedures  dis¬ 
cussed  above. 


10.  EXPERIMENTAL  RESULTS 

During  the  period  covered  by  this  report,  experi¬ 
mental  work  was  focused  mainly  on  the  following  areas. 

10.1  Quantitative  Absorption  Measurements  on 
Tri-n-Butyl  Phosphate"! 

Measurements  begun  in  the  previous  contract  period 
were  completed.  Initial  measurements  were  repeated  and  some 
additional  absorption  bands  were  investigated.  The  com¬ 
pleted  work  was  submitted  for  printing  in  the  original  re¬ 
port  (Survey  of  Infrared  Absorption  Phenomena  and  Experi¬ 
mental  Methods  for  Liquids  and  Solid  Films.  DAAK  10-77-C- 
0121,  Oct.  1977).  No  data  will  be  reported  here. 

10 . 2  Infrared  Reflection  Methods  for  Powders 

Developmental  work  was  continued  on  methods  of  deter¬ 
mining  the  optical  constants  using  reflection  methods.  Var¬ 
iable  angle  reflection  accessories,  obtained  for  this  study, 
were  mounted  in  the  spectrophotometer  for  test  measurements. 
Special  sample  holders  were  constructed.  Optical  alignment 
of  reflection  accessories  was  achieved.  It  was  found  that 

a  more  precise  angle  of  incidence  determination,  than  was 
available  from  the  stock  unit,  was  necessary.  A  He-Ne  laser 
was  purchased  and  used  for  the  angle  of  incidence  measure¬ 
ments.  Calibration  curves,  resulting  from  these  measure¬ 
ments,  were  obtained.  As  discussed  above  in  paragraph  9.1, 
preliminary  work  showed  that,  for  example,  carbon  and  gra¬ 
phite  powder  samples  could  be  cold  pressed,  using  a  die  and 
hydraulic  press,  to  yield  reasonably  good  solid  disks.  The 
surface  smoothness  of  the  pressed  disks  were  not  signifi¬ 
cantly  different  from  that  of  polished  solid  specimens, 
where  available.  Refinements  in  the  sample  holder  design 
were  found  to  be  necessary,  especially  for  the  carbon  disks 
which  tend  to  crumble  quite  easily.  Test  procedures  were 
begun  using  powders  and  solid  samples  of  known  optical 
constants. 


A  computer  interfaced  data  processing  system  was  also 
introduced  during  this  time.  The  system  is  showh  schemati¬ 
cally  in  the  block  diagram  below: 


Radiation  intensity  vs.  wavenumber  data  from  the  infrared 
spectrophotometer  is  digitized  by  a  digital  display  acces¬ 
sory.  An  interface  unit  encodes  the  data  and  directs  the 
Teletype  to  punch  the  coded  data  onto  paper  tape  while 
printing  out  the  intensity  values  to  four  significant  fig¬ 
ures  on  chart  paper.  The  paper  tape  is  then  transported  to 
the  computer  center  where  processing  and/or  calculations 
are  performed.  After  considerable  initial  difficulty  with 
the  analog  to  BCD  processing  components,  the  system  has 
been  rendered  functional.  A  primary  application  of  the 
system  is  to  extract,  by  computer  calculation,  the  values 
of  n  and  k  from  the  input  values  of  the  reflected  intensi¬ 
ties  and  angle  of  incidence.  The  inability  to  invert,  say 
I,  6  data  to  n,  k  data,  was  discussed  in  paragraph  6.1. 

See,  for  example,  Egns.  (58)  and  (59). 


35 


DISTRIBUTION  LIST  FOR  ARCSL-CR-80024 


Names 


Chemical  Systems  Laboratory 
Aberdeen  Proving  Ground,  MD  21010 

Office  of  the  Director 
Attn:  DRDAR-CLG 

CB  Detection  A.  Alarms  Division 
Attn:  DRDAR-CLC 

Developmental  Support  Division 
Attn:  DRDAR-CLJ-L 
Attn:  DRDAR-CLJ-R 
Attn:  DRDAR-CLJ-M 

Munitions  Division 
Attn:  DRDAR-CLN 
Attn:  DRDAR-CLN -S 

Physical  Protection  Division 
Attn:  DRDAR-CLW-P 

Research  Division 

Attn:  DRDAR-CLB 

Attn:  DRDAR-CLB-B 

Attn:  DRDAR-CLB-P 

Attn.  DRDAR-CLB-T 

Attn:  DRDAR-CLB-PS  (Mr.  Vervier) 

Attn:  DRDAR-CLB-PS  (Dr.  Stuebing) 

Attn:  DRDAR-CLB-PS  (Mr.  Frickel) 

Systems  Development  Division 
Attn.  DRDAR-CLY-A 
Attn:  DRDAR-CLY-R 

DEPARTMENT  OF  DEFENSE 

Administrator 

Defense  Technical  Information  Center 
ATTN:  DTIC-DDA-2 
Cameron  Station,  Bklg  5 
Alexandria,  VA  22314 


Office  of  the  Director 
Defense  Research  and  Engineering 
Attn:  Dr.  T.C.  Walsh,  Rm  3D-1079 
Washington,  DC  20310 

Institute  for  Defense  Analysis 
400  Army-Navy  Drive 
Attn:  L.  Biberman 
Attn:  R.  E.  Roberta 
Arlington,  VA  22202 

Advanced  Research  Projects  Agency 
1400  Wilson  Boulevard 
Arlington,  VA  22209 


Copies 


1 


3 

2 


1 

1 


1 

1 

1 

1 

10 

1 

I 

6 


12 


1 


1 

I 

1 


Names  Copies 

DEPARTMENT  OF  THE  ARMY 


Commander 

US  Army  Research  Office  Durham  I 

Box  CM,  Duke  Station 
Durham,  NC  27706 

HQDA  (DAMO-SSC)  I 

HQDA  (DAMA-ARZ,  Dr.  Verdcrame)  1 

HQDA  (DAMA-CSM-CM)  I 

HQDA  (DAMI-FIT)  I 

WASH  DC  20310 


US  ARMY  MATERIEL  DEVELOPMENT  AND 
READINESS  COMMAND 

Commander 

US  Army  Materiel  Development  and  Readiness  Command 
Attn:  DRCDE-DM 
Attn:  DRCLDC 
Attn:  DRCMT 
Attn:  DRCSF-S 
Attn.  DRCDL  (Mr.  N.  Klein) 

Attn:  DRCB1  (COL  Gearin) 

Attn:  DRCDMD-ST  (Mr.  T.  Shirata) 

S001  Eisenhower  Ave. 

Alexandria,  VA  22333 

Commander 

US  Army  Foreign  Science  &  Technology  Center 

Attn:  DRXST-MT-2  I 

Attn:  DRXST-CE  (Mr.  V.  Raguc)  I 

220  Seventh  St.,  NE 

Charlottesville,  VA  22901 


Commander 

US  Army  Missile  Command 
Redstone  Scientific  Information  Center 

Attn:  Chief,  Documents  1 

Attn:  DRDMI-CGA  (Dr.  B.  Fowler)  1 

Attn:  DRDMI-TE  (Mr.  H.  Anderson)  1 

Attn:  DRDMI-KL  (Dr.  W.  Wharton)  1 

Redstone  Arsenal,  AL  3S809 


US  ARMY  ARMAMENT  RESEARCH  AND 
DEVELOPMENT  COMMAND 

Commander 

US  Army  Armament  Research  and  Development  Command 
Attn:  DRDAR-TSS  S 

Dover,  NJ  07801 

US  ARMY  ARMAMENT  MATERIEL  READINESS  COMMAND 
Commander 

US  Army  Armament  Materiel  Readiness  Command 

Attn:  DRSAR-ASN  I 

Attn:  DRSAR-PE  1 

Rock  Island,  IL  61299 


37 


DISTRIBUTION  LIST  FOR  ARCSL-CR-80024  (Contd) 


Names  Copies 

Commander 

Harry  Diamond  Laboratories 

Attn:  DRXDO-RDC  (Mr.  D.  Giglio)  1 

2800  Powder  Mill  Road 
Adeiphi.  MD  20783 


Chief,  Office  of  Missile  Electronic  Warfare 
US  Army  Electronic  Warfare  Laboratory 

Attn:  DRSEL-WLM-SE  (Mr.  K.  Larson)  1 

White  Sands  Missile  Range,  NM  88002 

Project  Manager  for  Smoke/Obscurants 

Attn:  DRCP-SMK  2 

Aberdeen  Proving  Ground,  MD  21005 

Commander 

Atmospheric  Sciences  Laboratory 

Attn:  DRSEL-BR-AS-P  1 

Attn:  DRSEL-BR-MS-A  (Dr.  R.  Gomez)  1 

Attn:  DRSEL-BL-AS-DP  (Mr.  1.  Lindberg)  1 

Attn:  DRSEL-BL-SY  (Mr.  F.  Horning)  1 

While  Sands  Missile  Range,  NM  88002 

Director 

US  Army  Materiel  Systems  Analysis  Activity 

Attn:  DRXSY-D  (Dr.  Fallin)  1 

Attn:  DRXSY-MP  1 

Aberdeen  Proving  Ground,  MD  21005 

Director 

Night  Vision  Laboratories 

Attn:  DRSEL-NV-VI  (Mr.  R.  Moulton)  1 

Attn:  DRSEL-NV-VI  (Mr.  R.  Bergemann)  1 

Fort  Belvoir,  VA  23651 

US  Army  Mobility  Equipment  Research  and 
Development  Center 

Attn.  C  ode/DROME-RT  (Mr.  O.  F.  Kezer)  1 

I  ort  Belvoir.  VA  22060 

Commander 

US  Army  Electronics  Command 

Attn  DRSEL-CT-LG  (Dr.  R.  G.  Rohde)  1 

Aitn:  DRSEL-CT-l  (Dr.  R.  G.  Buser)  1 

Attn:  DRSEL-WL-S  (Mr.  J.  Charlton)  1 

Fort  Monmouth,  NJ  07703 

Director 

US  Army  Ballistic  Research  Laboratories 

Attn:  DRXBR-DL  (Mr.  T.  Finnerty)  1 

Attn:  DRXBR-P  (Mr.  N.  Gerri)  1 

Attn:  DRDAR-BLB  (Mr.  R.  Reitz)  1 

Attn:  DRDAR-BLB  (Mr.  A.  LaGrange)  1 

Aberdeen  Proving  Ground,  MD  21005 

Cdr.  APG 

ATTN:  STEAP-TL  I 

Aberdeen  Proving  Ground,  MD  21005 


Names  Copies 

Commander 

US  Army  Test  t  Evaluation  Command 

Attn:  DRSTE-FA  1 

Aberdeen  Proving  Ground,  MD  21005 

Commander 

Dugway  Proving  Ground 

Attn:  STEDP-PO  1 

Attn:  Technical  Library,  Docu  Sec  I 

Attn:  STEDP-MT-DA-E  I 

Attn:  STEDP-MT  (Dr.  L.  Salamon)  I 

Dugway,  UT  84022 


US  ARMY  TRAINING  &  DOCTRINE  COMMAND 

Commandant 

US  Army  Infantry  School 

Combat  Support  &  Maintenance  Dept 

Attn:  NBC  Division  I 

Fort  Benning,  GA  31905 

Commandant 

US  Army  Missile  &  Munitions  Center  &  School 

Attn:  ATSK-CD-MD  I 

Attn:  ATSK-DT -MU-EOD  1 

Redstone  Arsenal,  AL  35809 

Commander 

US  Army  Logistics  Center 

Attn:  ATCL-MM  1 

Fort  Lee,  VA  23801 

Commander 
HQ,  USA  TRADOC 

Attn:  ATCD-TEC  (Dr.  M.  Pastel)  I 

Fort  Monroe,  VA  23651 

Commander 

US  Army  Ordnance  &  Chemical  Center  &  School 

Attn:  ATSL-CD-MS  1 

Aberdeen  Proving  Ground,  MD  21005 


DEPARTMENT  OF  THE  NAVY 
Commander 

Naval  Surface  Weapons  Center 

Attn:  Tech  Lib  &  Info  Svcs  Br  I 

White  Oak  Laboratory 
Silver  Spring,  MD  20910 

Commander 

Naval  Intelligence  Support  Center  1 

4301  Suitland  Road 
Washington,  DC  20390 

Commander 

Naval  Surface  Weapons  Center 
Dahlgren  Laboratory 

Attn:  DX-21  I 

Dahlgren.  VA  22448 


DISTRIBUTION  LIST  FOR  ARCSLCR-80024  (Contd) 


Name* 

Copies 

Names 

Commander 

DEPARTMENT  OF  THE  AIR  FORCE 

Naval  Weapon*  Canter 

Attn:  Code  3311  (Dr.  R.  Bird) 

1 

HQ,  Foreign  Technology  Division  (AFSC) 

Attn.  Code  382  (Dr.  P.  St,  Amand) 

1 

Attn:  PDRR 

Attn:  Code  3122  (Dr.  Hindman) 

1 

Wright-Patterson  AFB,  OH  45433 

China  Lake.  CA  935SS 

Commanding  Officer 

Commander 

Armament  Development  4  Test  Center 

Naval  Weapon*  Support  Center 

1 

Attn:  DLOSL  (Technical  Library) 

Attn:  Code  5041  (Mr.  D.  Johnaon) 

Eglin  AFB,  FL  32542 

Attn:  Code  5042  (Mr.  C.  Dtaerman) 

1 

Crane.  IN  47522 

Commander 

Commander 

tKFiimM 

Attn:  ECJS-O/LTC  James  H.Alley 

Naval  Research  Laboratory 

1 

APO,  NY  09128 

Attn:  Code  5709  (Mr.  W.  E.  Howell) 

4555  Overlook  Avenue,  SW 

USA  BIOMEDICAL  LABORATORY 

Washington,  DC  20375 

Attn.  SGRD-UV-Za 

US  Army  Trailing  A  Doctrine  Command  Liaison 

German  NBC  School 

clo  7261  MUNSS,  Box  2« 

APO  New  York  09035 


